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Electron-positron bremsstrahlung and pair creation in very 
high magnetic fields 



o 

(N 

wo: 
< 

X 

6 

in 

c3 



P. B. Jones* 

Department of Physics, University of Oxford, Denys Wilkinson Building, 
Keble Road, Oxford OX1 3RH 



> 

m 
co 

q 

00 

o 
o 



X 



ABSTRACT 

Cross-sections for Rutherford scattering, Coulomb bremsstrahlung and pair creation, 
have been calculated at very high magnetic fields in order to investigate the photo- 
production of protons at the polar caps of pulsars whose spin is antiparallel with 
the polar magnetic flux density. The Landau-Pomeranchuk-Migdal effect at very high 
magnetic fields is included in a simple electron Green function. 
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1 INTRODUCTION 

Work on quantum electrodynamic processes at magnetic flux 
densities of the order of the critical field B c = 



4.41 x 10 13 G has been reviewed recently, and very com- 
pletely, by Harding & Lai (2006). The processes concerned 
are primarily those relevant to collisions within plasma at 
low altitudes above the surface of a neutron star and there- 
fore do not include relativistic Coulomb bremsstrahlung and 
pair creation in the interaction of electrons and photons 
with ions. There appear to be no published cross-sections 
for these processes at B ~ B c apart from a number of pa- 
pers on the bremsstrahlung energy-loss of non-relativistic 
electrons. But of these, only Lieu (1981) and Lauer et al 
(1983) have treated bremsstrahlung as a process of second 
order in the electron-photon coupling and the work of the 
latter authors is also limited to collisions in which both ini- 
tial and final electrons are in the lowest Landau state. Bus- 
sard (1980) and Langer (1981) have considered photon pro- 
duction in non-relativistic electron-ion plasma collisions but 
have treated it as a sequence of first-order processes, that 
is, Coulomb excitation of a higher Landau state followed by 
cyclotron emission. 

However, the reverse flux of electrons or positrons pro- 
duced by pair creation in the open magnetic flux-line re- 
gion above neutron-star polar caps must produce electro- 
magnetic showers in the condensed matter at the surface 
whose properties are determined by the relativistic second- 
order processes. In neutron stars with spin fl such that 
tt ■ B > 0, the outward particle flux consists of electrons so 
that nuclear reactions within the positron-initiated showers 
have no obvious observable effect. But in the opposite case, 
fi-B < 0, the outward particle flux is positively charged and, 
with positrons, also includes protons that are produced in 
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the shower and diffuse to the surface. It is likely that the 
temporal behaviour of these different components does lead 
to instabilities that are the basis of observable phenomena 
(Jones 2010). In general, there appears to be no reason a 
priori why both ft ■ B cases should not be present in the 
neutron-star population and it is of interest to understand 
how each might be observed in the electromagnetic spec- 
trum. Estimates of shower development in previous work 
relied on the correspondence principle, that for large values 
of the Landau quantum number the zero-field cross-sections 
are valid, but this is unlikely to be satisfactory at B ~ B c 
where the number of contributing Landau states can be 
small. The present paper attempts to give cross-sections that 
are rather better than order-of-magnitude estimates but are 
not of high accuracy owing to the discrete nature of the 
Landau spectrum and to truncations that are made in per- 
forming the calculations. Its purpose is to see whether and 
in what way the cross-sections differ qualitatively from those 
at zero field. Nucleon production in electromagnetic showers 
occurs almost entirely through formation and decay of the 
giant dipole resonance. Thus the paper is concerned, par- 
ticularly, with factors that influence the low-energy photon 
track length and its depth distribution. 

With neglect of radiative corrections and of the anoma- 
lous magnetic moment, the energy states of an electron in a 
uniform magnetic field are 



E = ±y/p 2 + l + 2nB 



(1) 



in which p is the longitudinal momentum component, paral- 
lel with B, and n = 0, 1, 2, ... is the Landau quantum num- 
ber. Apart from the Appendix and where explicitly stated 
otherwise, momentum and energy are here expressed in units 
of mc and mc 2 ; the magnetic flux density B is in units of 
B c . In electron- photon interactions, the conserved quanti- 
ties are energy and the longitudinal momentum component. 
Thus the kinematic behaviour of an electron is that of a par- 
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tide confined to one dimension and of mass m^Jl + p 2 ± in 
which we have defined, purely for ease of notation, a notional 
transverse momentum component p± — V2nB. 

We have adopted the solutions of the Dirac equation 
found by Johnson & Lippmann (1949) and give a summary 
of calculational details in the Appendix. These are unre- 
markable except that owing to the one-dimensional kine- 
matics of electrons and the small number of Landau states 
that contribute significantly at B ~ B c , we have found it 
convenient, given that numerical computation must be in- 
volved, not to proceed with the usual approach of covariant 
formalism and Green functions but instead to write down 
the transition matrix to second order directly in terms of 
time-ordered matrix elements calculated using an explicit 
representation of the Dirac matrices. It is also appropriate 
to mention here that the states given by equation (1) are 
two-fold degenerate for n > only if radiative level shifts, 
natural widths, and the anomalous magnetic moment of the 
electron are neglected. But for the high-energy spin-averaged 
or summed processes considered here, the fine structure is 
of no consequence and we use only the completeness of the 
Johnson-Lippmann eigenfunctions. References and some fur- 
ther details are given in the Appendix. 

The most simple problem here is that of finding the re- 
placement for Rutherford's scattering formula. Small-angle 
Coulomb scattering in the zero-field case is replaced at 
B ~ B c by the Coulomb excitation of Landau states. This 
has been treated by Bussard (1980) and Langer (1981) 
but only in the non-relativistic limit. The relativistic cross- 
sections that we require appear not to have been published 
previously and are given in Section 2. There have been many 
calculations of transition rates for cyclotron emission, and 
a summary is contained in the review of Harding & Lai 
(2006). New transition rates are given in Section 3 to al- 
low comparison with previous work and, more specifically, 
to obtain the distribution of photon transverse momentum 
k± from the decay of n > 1 Landau states which is of 
crucial importance in the development of electromagnetic 
showers at B ~ B c . Cross-sections for Coulomb pair cre- 
ation and bremsstrahlung are obtained in Sections 4 and 5. 
In the case of pair creation, there appears to be no previ- 
ously published work. For bremsstrahlung, the work of Lieu 
(1981) and Lauer et al (1983) is in the non-relativistic limit 
so that comparison is not possible. There is also the fur- 
ther problem that calculations relying on truncation of the 
number of Landau states lack external tests of their correct- 
ness because the zero-field limit is not available. But where 
minor tests of correctness exist, we find that they are sat- 
isfied. Finally, with primary electron energies of the order 
of 10 GeV and the high condensed matter densities pre- 
dicted at the neutron star surface, it is essential to consider 
the Landau-Pomeranchuk-Migdal effect and its influence on 
Coulomb pair creation and bremsstrahlung cross-sections. 
Its existence and broad properties at zero field have been 
verified through experiments at particle accelerators. We re- 
fer, in particular, to the recent work of Hansen et al (2004) 
and to the review of Klein (1999). We find that the effect 
exists at B ~ B c and that the one-dimensional electron kine- 
matics makes possible a particularly simple treatment which 
is given in Section 6. The cross-section calculations of Sec- 
tions 4 & 5 have been repeated with allowance for the LPM 
effect which is shown to significantly change early shower 



development. The final section gives a qualitative summary 
of the cross-sections obtained and of the effect of magnetic 
fields of the order of the critical field on shower development. 



2 RUTHERFORD SCATTERING 

Small-angle Coulomb scattering at zero field, in which the 
transverse momentum component is a continuous variable, 
is replaced by a Landau transition n — > n' with longitudinal 
moment transfer q. Free atoms in high magnetic fields have 
an axially symmetric but complex electron density distribu- 
tion. But the departures from uniform electron density in 
condensed matter are not large except for a small number 
of inner orbitals and we therefore adopt the potential ZeV 
of a point charge within a spherical Wigner-Seitz unit cell 
of uniform electron density. Its radius is 

r ws = 2.15 x io- 10 Z°- 23 B- - 40 cm, (2) 

obtained from the ion number density calculated by Medin 
& Lai (2006) for atomic number Z. Calculational details arc 
summarized in the Appendix. The symmetry of V indicates 
that Landau functions in cylindrical polar coordinates are 
the appropriate choice. The quantum number I denotes the 
spatial degeneracy of these states and the procedure is to 
average the transition rate over all I in the range < / < Z m . 
Each state has a guiding centre radius y/2l + Irs, where tb 
is the cyclotron radius. The incident flux is then c divided by 
(2l m + l)nr B , and the cross-section at incident momentum 
p >• 1, with no spin- flip, is the limit l m — > oo of, 




J2\{p + q,n',l'\V\p,n,l)\ 2 . (3) 

1=0 

This expression is greatly simplified by the selection rule 
5(n — I) — arising from equation (A3) and from the ax- 
ial symmetry of V. Spin-flip cross-sections are smaller by a 
factor of the order of (p' ± — p±) 2 /4:E 2 relative to unity and 
as our interest is primarily in electrons with p 3> mc they 
have not been calculated. The range of initial-state quan- 
tum numbers I should be so large that the cross-section is 
independent of l m in the case of an isolated atom. We have 
confirmed that, for B — B c , this condition is satisfied ad- 
equately for l m « 30 at which value guiding centres are 
at radii approaching the Wigner-Seitz cell radius. Numeri- 
cal cross-sections a nn i are given in Table 1, for n, n' < 4. 
They are in units of lbn = 10~ 24 cm 2 and are for Z = 26 
and longitudinal momentum p = 20, but are almost exactly 
independent of p. Calculation to the lowest order in V is 
satisfactory here, as for the zero-field small-angle Ruther- 
ford formula, owing to the pole in the amplitude at zero 
momentum transfer. Thus the cross-sections are linear in 
Z 2 . 

The cross-sections of Table 1 are valid in the relativistic 
limit of p ^> 1 so that comparison with the near-threshold 
cross-sections found by Bussard (1980) and Langer (1981) 
is not possible. But the standard Rutherford cross-section 
at small angles can be expressed in terms of a continu- 
ous variable n through the correspondence principle relation 
p± = \/2nB. It is dan/dn = 0.50Z 2 /n 2 B bn. The numerical 
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Table 1. Cross-sections in units of 10~ 24 cm 2 for Coulomb- 
induced Landau transitions n — > n' , with n, n' ^ 4, arc almost 
exactly independent of longitudinal momentum at p ^> mc but 
have been calculated for electrons of p = 20 mc incident on an 
isolated neutral atom of atomic number Z = 26. The magnetic 
flux density in the first column is in units of B c = 4.41 X 10 13 
G. The significance of the n — ¥ n transitions is considered in the 
text and later in Section 6. 
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values and dependences on n and B are quite close to those 
given in the Table for an initial state with n — 0. The n — > n 
cross-sections are well-defined by equation (3) even though 
they are to an evanescent state whose asymptotic final 
form with q = is indistinguishable from the unscattered. 
These Sn = transitions to virtual states of finite q ap- 
pear in the calculation of second-order processes, in partic- 
ular bremsstrahlung and pair creation cross-sections, and in 
condensed matter are relevant to the Landau-Pomeranchuk- 
Migdal effect considered in Section 6. 

The cross-sections in Table 1 are typically two orders 
of magnitude larger than those we obtain for Coulomb 
bremsstrahlung in Section 5. Thus Rutherford scattering 
excites higher Landau states whose decay is an important 
source of photons, particularly in the later stages of shower 
development. We shall refer to this in greater detail in Sec- 
tion 7. 



k sin 2 9 = E - (El - sin 2 8{pl - p't)) 1/2 , 



(4) 



in which Eq = U 1 + p\ . The partial transition rate for 
unpolarized electrons is 



E' 



1 f 1 k 2 

r„„' = -r / d(cos6) — — — — 

n j j 2irh 6 c E' + ck cos 2 8 



^j-jT-y^y^ \{l',n',s',\i,e\ea.-k\n,l,s)\ 2 , (5) 



2(l m + 1) 



e s,s' l,l'=0 



summed over photon polarization states eii and e±, respec- 
tively parallel with and perpendicular to k x B, and over 
final electron spins. The radiation field vector potential is A 
given by equation (A6) and the final-state electron energy is 
E' . The transition rate is expressed here as an average over 
lm + 1 spatially degenerate initial states and we have ver- 
ified that the computed rates are Z m -independent, as they 
should be, and have used the values l m = and l' m — 12. 
The computed values, in units of W~ 3 uib, where lob is the 
cyclotron angular frequency, are given in Table 2 for n ^ 4. 
The summed transition rates given in the final row of each 
section are the sole results of this paper that can be com- 
pared directly with previously published work. Agreement 
with the rates shown in Fig. 3 of Harding & Lai (2006) is 
satisfactory insofar as accurate comparison is possible. It is 
interesting to see that there is effectively no selection rule 
on Sn and that the partial rates display the <5n-dependence 
noted by these authors in that, with increasing magnetic flux 
density, there appears an enhanced rate for direct transitions 
to the n — ground state. The dependence on photon po- 
larization is not given here. It is unremarkable: transitions 
to ey are approximately twice as strong as those to the e±. 
Photon angular distributions vary with Sn but are also un- 
remarkable and are not given here. For Sn — —4 they have 
a maximum at 6 = 7r/2; for Sn = —1 the maxima are near 

e = o,7T. 

The maximum photon transverse momentum occurs at 
6 = 7r/2 and is, 



k±. 



Pi 



(6) 



3 CYCLOTRON EMISSION 

Transition rates for relativistic cyclotron emission are obvi- 
ously relevant to the low-altitude plasma above polar caps 
at B ~ B c and therefore have been calculated by a number 
of authors (Herold, Ruder & Wunner 1982, Latal 1986, Bar- 
ing, Gonthier & Harding 2005; see also the review of Harding 
& Lai 2006). The principal reason for our calculation here, 
which is limited to unpolarized electrons, is to investigate 
the sequence of partial transition rates by which an electron 
in Landau state n decays to n = 0. These determine the dis- 
tribution of the transverse photon momentum k± — ks'm9, 
where 8 is the photon angle with respect to B, and hence 
electromagnetic shower development at B ~ B c . 

In the rest frame of the initial electron, the relation 



It is not possible to make a compact statement about the dis- 
tributions of k± because the values of k± imax are functions of 
both n and B. The critical thresholds for single-photon mag- 
netic pair creation are those at k± — 2 for electron-positron 
Landau numbers n± — n T = and k± = 1 + Vl + 2B 
for n± — with n T = 1. Rates immediately above the lower 
threshold are partially inhibited by a selection rule but above 
the higher threshold, transition rates for magnetic pair cre- 
ation are so high that photon mean free paths are small 
compared with the radiation length, which is the basic unit 
of length for shower development. But Table 2 shows that, 
except for the highest magnetic field, the majority of tran- 
sitions from n ^ 4 have k± < 2 and so are not followed by 
magnetic pair creation. Thus cyclotron emission, following 
excitation by Rutherford scattering, is an important source 
of photon track length in a shower. 
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Table 2. Rates for cyclotron emission with Landau transitions 
rt — > n' are given in units of 1.0 X 10 _3 ws for unpolarized elec- 
trons. The sums of the partial transition rates for n < 4 are in 
the final row of each section. 
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sum 


0.94 


1.41 
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0.35 
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0.39 


0.22 
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0.46 


0.25 


0.17 







0.62 


0.41 


0.34 


0.29 




sum 


0.62 


0.87 


0.98 


1.03 



4 PAIR CREATION 

Coulomb pair creation in a uniform magnetic field requires 
a longitudinal momentum transfer q from the nucleus, 

q = k(l- cos6) - J- (1 +pi+) - ^- (1+PX-) , (7) 

an expression valid for final-state electron and positron mo- 
menta p± ^> 1, and therefore differs from the zero-field case 
in being a function of both k and k± . The momentum trans- 
fer is a minimum for p + — p- and p±+ = p±- — for which 
conditions it falls to zero at k± = 2, the kinematic threshold 
for magnetic pair production. The transition matrix element 
to second order for an isolated atom is, 

M[± = (p-,n-,l-,s-\ - Ze 2 V\p- - q,n-,l-,s-) 

(E-(p-) - E-(p- -q)+ir,r 1 

(P- - q,n-,l-,s-\ea • A| -p + ,n+,l+, -s+) 

+ (P+, n+, l+, s+\Ze 2 V\p+ - q, n+, l+, s+) 

(E + (p + )-E+(p + -q) + iri)- 1 

(p-,n-,l-,s-\ea ■ A| - p+ + q,n + ,l+,-s+) (8) 

in which rj > is infinitesimal. The subscripts || and _L refer 
to the photon polarization with respect to k x B. The spin 
and spatial degeneracy states are labelled by s± and l± (see 
equations A7 and A8). There are two further time-ordered 
terms in the amplitude but they involve virtual pair creation 
in the Coulomb field. They have energy denominators that 
are large, of the order of k rather than k~ as in the above, 
and for k 1 they can be neglected. The total cross-section 
is then, 

a P ± = Z 2 (^-Y(A-) 2 f-y/^dE- E-E + 
l± \hcj Kmc) J ^/ 1+p 2 ± _ 2tt p-p + 



EEEK-1 2 <°) 

for pair creation by a photon of momentum k, transverse 
momentum k±, and of polarization ey or ej_. In this expres- 
sion, the charges have been factored out from the matrix 
element (8). There are a number of approximations here 
that merit some consideration. Firstly, truncation of the 
Landau quantum numbers to n± ^ 4 has been tested by 
examination of the computed cross-section as a function of 
n m for < n m ^ 4. The basis for this truncation is that, 
for fixed k, the amplitude given by equation (8) is approx- 
imately oc pX+, also that matrix elements tend rapidly to 
zero as they are constructed between states with increas- 
ingly different quantum numbers and hence nodal surfaces. 
There is evidence of reasonably satisfactory convergence, but 
cross-sections obtained here should be seen as lower limits. 
Truncation of the sum over the spatial degeneracy quantum 
numbers to l± < 12 again produces fair convergence. But in 
this case, there is a natural upper limit determined by the 
radius of the Wigner-Seitz cell. Finally, the matrix elements 
of the Coulomb potential ZeV have been limited to those 
with 8n = 0. This is satisfactory, for the purposes of the 
present work, because these are typically much larger than 
off-diagonal matrix elements, as shown by the cross-sections 
of Table 1. But errors arising from this latter truncation are 
not positive-definite. 

Numerical values of the isolated-atom total cross- 
sections are given in Table 3 for ranges of values of B, k and 
k± and for both photon polarizations. They are for Z — 26 
but being of second order in electron-photon coupling, are 
linear in Z 2 . The differential cross-sections are unremark- 
able symmetric functions of the final-state energies E± , not 
strongly peaked within the kinematic limits present in the 
integral of equation (9). They are not given here primar- 
ily because in the relevant case of condensed matter they 
may be unpredictably distorted. The explanation for this is 
a consequence of the partial order which may be present in 
neutron-star surface condensed matter, and is deferred until 
Section 7. 

A problem in this Section is that there are no external 
tests of the validity of expressions derived from equations (8) 
and (9) or of their numerical evaluation. But we note that 
the differential cross-sections are symmetric in E± and that 
under the restriction n m = 0, the cross-section is cry 3 = 
for photons polarized parallel with k x B. This is consistent 
with the fact that for this polarization, the matrix element 
(A7) should be exactly zero for n+ = n_ = (see Semionova 
& Leahy 2001). 

The cross-sections in Table 3 are quite slowly varying 
functions of k and can be compared with the asymptotic 
k ^> 1 Bethe-Heitler value of a F — 5.0 bn for unpolar- 
ized photons incident on a screened Z = 26 nucleus in zero 
magnetic field. The correspondence principle indicates that 
the Bethe-Heitler expression with modified screening should 
be a fair approximation at moderate fields, of the order of 
10 12 G, for which the Landau quantum number associated 
with p± ~ 1 is large. It is difficult to obtain cross-sections 
for fields between 10 12 G and those of the Table because 
the truncations we have used are unlikely to be satisfac- 
tory approximations. But as a function of decreasing B, the 



5 



Table 3. Numerical values of the total cross-section (in units 
of lbn = 10~ 24 cm 2 ) for pair creation in the Coulomb field of 
an isolated atom with nuclear charge Z = 26 are given for both 
photon polarizations in columns 3-5 of the upper sector of the 
Table. The photon momentum k and transverse momentum k± 
are in units of vac. The magnetic flux density is in units of the 
critical field B c = 4.41 X 10 G. The three-row lower sector gives 
cross-sections re-calculatcd in condensed matter with inclusion 
of the Landau-Pomcranchuk-Migdal effect. We refer to Section 6 
for the details of this latter calculation, which assumes a static 
structure function S(q) = 1. 
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0.24, 0.24 
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0.57, 1.11 




4000 


0.25, 0.25 


0.38, 0.41 


0.59, 1.13 
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40 


0.06, 0.06 


0.07, 0.08 


0.08, 0.23 




400 


0.07, 0.07 


0.08, 0.10 


0.11, 0.35 




4000 


0.07, 0.07 


0.09, 0.11 


0.11, 0.36 
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40 


0.61, 0.59 


1.22, 1.23 


1.99, 2.47 




4000 


0.21, 0.21 


0.51, 0.50 


0.78, 0.79 




40000 


0.02, 0.02 


0.05, 0.05 


0.08, 0.08 



cross-sections in the Table are not inconsistent with smooth 
convergence at intermediate fields toward the Bethe-Heitler 
value. The most noticeable features of the Table are: (i) 
the rapid increase of cross-section as k± approaches the n±- 
dependent thresholds for magnetic pair production at which 
q and hence the energy denominators in equation (8) be- 
come zero; and (ii), the rapid decrease of cross-sections at 
B > B c . 



5 BREMSSTRAHLUNG 

The emission of a photon of momentum k at an angle 9 
with B, and a change of Landau quantum number n — > n , 
requires a longitudinal momentum transfer from the atom 
of 

q = _ k{1 _ cose) + l+^_l+ll j (10) 

where the initial electron longitudinal momentum is p ^> mc 
and the final momentum is p' — p + q — k cos 6. 

The transition matrix element, to second order for an 
isolated atom, is 

M^ ± = (p , k, n , l' , s'\eoc ■ A\p + q, n, I, s) 
(E(p)- E(p + q) + i V )- 1 
{p + q,n,l,s\ - Ze 2 V\p,n,l,s) + 
(p , n , l',s'\ — Ze 2 V\p — q,n',l', s') 
(E(p') - E(p' - q) + ir,)- 1 
(p - q,k,n',l',s'\ea ■ A\p,n,l,s), (11) 



in which the subscripts || and _L again refer to the pho- 
ton polarization with respect to k x B. As in the previous 
Section, and for the same reason, we have neglected two 
time-ordered terms in the amplitude that involve virtual 
pair creation in the Coulomb field. The existence of single- 
photon magnetic pair creation thresholds implies that the 
bremsstrahlung cross-section in a high magnetic field is best 
expressed in terms of the independent variables k and k±. 
The cross-section defined by the energy-loss rate (equivalent 
to the radiation length) is then given by an integration over 
photon energy, 

° =J dk pHk-' < 12 > 

in which the differential bremsstrahlung cross section for 
an electron in Landau state n, summed over both photon 
polarization states, is, 




EEEEKI 2 - ( 13 ) 

l.V— e s,s' n' — 

The summations have been truncated to n m = 4 and 
l m = 12, as in the case of pair creation, but this approxi- 
mation is much less satisfactory for bremsstrahlung and has 
lead us to omit values for B < B c . Bearing in mind the 
importance of the variable k±, the cross-section cr Rad given 
by equations (12) and (13) has been divided into two com- 
ponents: af ad is for the interval < k± < 2.5 and a^ ad 
for 2.5 < k± < 5.0 mc. The computed values in units of 
1CP 24 cm 2 , for electrons in an initial n = Landau state 
incident on isolated atoms are given in the upper sector of 
Table 4. They are summed over both photon polarizations, 
but we note that the e± cross-section is typically an order 
of magnitude larger than that for ey . The computed values 
are almost completely independent of p and there is little 
cross-section outside the intervals of k± considered. The to- 
tal cross-section is simply denned here as the sum of trf^- 
Cross-sections for initial Landau states n > are of similar 
magnitude but are not given here. These initial states are, 
of course, naturally unstable against photon decay, which 
appears as a zero in the denominators of equation (11) and 
is discussed further in Section 7. 



6 THE LANDAU-POMERANCHUK-MIGDAL 
EFFECT 

In both this and previous work, we have adopted expressions 
for the zero-pressure ion number density, 

N = 2.5 x lO 28 ^ ' 7 ^ 1 ' 2 cm~ 3 , (14) 

and for the interatomic separation in a linear chain, parallel 
with B, 

a s « 1.5 x W- 10 Z 1/2 B- 1/2 cm, (15) 

conveniently summarizing the calculated values of Medin & 
Lai (2006). We anticipate that the bremsstrahlung and pair 
creation cross-sections obtained in previous Sections for iso- 
lated atoms will be considerably modified in the condensed 
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Table 4. Numerical values of the cross-section equivalent to the 
radiation length, defined by equation (12), are given in the up- 
per sector of the Table for electrons in the n = Landau state 
incident on isolated atoms of Z = 26. The values of electron mo- 
mentum p are in units of mc, and B is in units of the critical 
field B c = 4.41 X 10 13 G. Columns 3 and 4 are partial cross- 
sections defined for the photon transverse momentum intervals 
< kj_ < 2.5 and 2.5 < k± < 5.0 mc, respectively. The total 
cross-section in column 5 is defined here as the sum of these par- 
tial cross-sections. The 15 rows in the lower sector of the table 
give the above cross-sections re-calculated in condensed matter 
assuming a static structure function S(q) = 1 but with inclusion 
of the Landau-Pomcranchuk-Migdal effect. We refer to Section 6 
for further details. 



p 


B 


~Rad 
°1 


—Bad 
°2 


g-Rad 


200 


1.0 


4.83 


1.85 


6.68 




3.0 


0.77 


1.44 


2.21 




10.0 


0.08 


0.25 


0.33 


2000 


1.0 


4.82 


1.85 


6.67 




3.0 


0.77 


1.44 


2.21 




10.0 


0.08 


0.25 


0.33 


20000 


1.0 


4.82 


1.85 


6.67 




3.0 


0.77 


1.44 


2.21 




10.0 


0.08 


0.25 


0.33 


200 


1.0 


22.83 


1.87 


24.70 




3.0 


9.94 


1.63 


11.54 




10.0 


3.92 


0.76 


4.68 


2000 


1.0 


20.04 


2.41 


22.45 




3.0 


6.95 


3.59 


10.54 




10.0 


1.75 


1.88 


3.63 


20000 


1.0 


5.35 


2.15 


7.50 




3.0 


1.32 


2.66 


3.98 




10.0 


0.22 


0.78 


1.00 


60000 


1.0 


1.52 


0.95 


2.47 




3.0 


0.32 


1.05 


1.37 




10.0 


0.05 


0.25 


0.30 


200000 


1.0 


0.23 


0.21 


0.44 




3.0 


0.04 


0.21 


0.25 




10.0 


0.01 


0.04 


0.05 



matter of the neutron-star surface. In particular, it is essen- 
tial to consider the Landau-Pomeranchuk-Migdal (LPM) ef- 
fect, whose existence at zero field has been established with- 
out doubt by experiments with laboratory electron beams of 
several hundred GeV (Hansen et al 2004). The essence of the 
effect is that, for example in bremsstrahlung and at high en- 
ergies, the isolated-atom longitudinal coherence length (of 
the order of g _1 ) can be so long that the coherence that 
would be present in the case of an isolated atom is par- 
tially removed by multiple scattering of the initial or final- 
state electron within the medium. We refer to the review by 
Klein (1999) for further details, and for a discussion of other 
processes in condensed matter, such as the dielectric effect, 
that are not relevant here. In the present Section, we exam- 
ine the LPM effect at high magnetic fields and find that the 
one-dimensional nature of kinematics at B ~ B c allows a 
particularly simple description. 

We begin by considering the effect of scattering of an 
electron with initial longitudinal momentum p ^> mc by a 
(dimensionless) potential U l = —Ze 2 V /he whose origin is 
at z = Zi. Wc consider transitions with no change of Landau 
quantum number so that the final state is given by the Green 



function, 

G = -i0(t - t') ^(ri)<(^)e lp ' (z - z '»- B ' (, - t '», 

p'>0 



(16) 



constructed from the Landau functions given by equation 
(A3). The truncation to Sn — is possible because, as wc 
found in Section 2, the Sn ^ transition rates are at least 
an order of magnitude smaller. It also excludes spin-flip and 
back-scattering (p' < 0) which are both negligible in the 
present context. The scattered state is 

J d z r'dt'G{r 1 t;r',t l )U i i> n {v'i_)e ipz '- im ' = 



f 

J — c 



dq Ug 
2ir —q + if) 



(17) 



in which U q is the matrix element of U % between states p and 
p' = p + q, and n > is infinitesimal. It is reached by intro- 
ducing the integral representation of the Heaviside function 
9, followed by integration over r' and t' . Equation (17) shows 
that the electron has components of momentum p + q and 
that these are significant at distances z — Zi < g -1 . Either by 
performing the integral, or taking the limit z — > oo, we can 
see that the prefactor outside the brackets evolves to an am- 
plitude Uq, as expected. Apart from the forward-scattered 
component of the wave- function given by equation (17), flux 
conservation shows that there is also an undisturbed compo- 
nent of amplitude (1 - \Uo\ 2 ) 1/2 . Equation (17), giving the 
forward-scattered state, is no more than a concrete expres- 
sion of the uncertainty principle, but it shows that forward 
scattering by a sequence of atoms, though unobservable in 
the motion of the electron itself, has the effect of broadening 
the distribution of q and so reducing the small-g amplitudes 
that mediate most of the bremsstrahlung cross-section at 
high energies. 

The simplicity of equation (17) enables us to find the 
Green function for an electron undergoing a sequence of 
forward scatters and so provides an easy way to calculate 
bremsstrahlung cross sections including the LPM effect. We 
begin by noting that the distribution in the variable q in 
equation (17) is unchanged in form by multiple forward scat- 
tering as is shown by the amplitude for a second and suc- 
cessive scattering at Zj, 



du I 
2^ / 

-OO J — C 



iq (z — zt) 



2tv —u + if] 



-Ul, 



dq 



q' + tn 

e iq(z-Zj) 



-q + if) 



(18) 



in which u = q — q' . We assume that the properties of U 
allow the limits of integration to be extended to ±oo with 
completion of the contour in the upper half of the complex-q' 
plane. We see that the amplitude has the spatial dependence 
z — zj of the last scattering. 

Consequently, the amplitude for an electron propagat- 
ing without back-scattering or Sn 7^ transitions in a long 
interval of z must be a superposition of g-distributions, each 
of the form of equation (17), with components which are 
necessarily identical except that their individual amplitudes 
decay exponentially as functions of z — Zj, reflecting a re- 
duction factor (1 — (l^oj 2 )) 1 ^ 2 for each atom traversed. This 
exponential representation is satisfactory provided the mean 
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square value of |?7q| 2 , satisfies the condition (|£/o| 2 ) <S 1, as 
it does here. Therefore, we can find the Green function that 
couples the electron to a final state of total longitudinal mo- 
mentum K« . It is given by, 



dze-^J_ 



k u z f°° dq U 3 q 



a ipz + iq(z — Zj ) — (z — Zj)/2\ 



2% —q + ir\ 

-1 



in which A = (a)/(|?7o| 2 }, where (a) is the mean interatomic 
separation Zj — Zj-i and we have integrated first over z and 
then over q by completion of the contour in the semi-infinite 
upper half of the complex q-plane. The right-hand side of 
equation (19) contains, within brackets, the momentum- 
space Green function for an electron propagating without 
back-scattering, spin-flip or Sn 7^ transitions. Forming the 
square modulus of the remaining term (and neglecting any 
inhomogeneity in the Uq ) yields the (one-dimensional) static 
structure function for a linear chain of N a atoms, 



3 = 1 



(20) 



which contains, as a function of q = Ku — p, the dependence 
of the transition rate on the condensed-matter structure. 
Allowance for transitions with change of Landau quantum 
number should not be difficult but is not essential because, 
as can be seen from the Table 1 cross-sections, their mean 
free path is considerably longer than A. 

Pair production and bremsstrahlung cross-sections in- 
cluding the LPM effect are obtained immediately, assuming 
a random distribution of atoms defined by a structure func- 
tion S(q) — 1, simply by substituting this Green function 
for the energy denominators contained in equations (8) or 
(11). We have done so and have repeated the calculations 
by which the Table 3 and 4 cross-sections were obtained. 
The revised cross-sections given in the lower sectors of these 
Tables are for Z — 26 and require some comment. 

In the case of pair production, the re-calculated cross- 
sections are given in Table 3 just for B = B c and show 
that, owing to the very high matter density (~ 2 x 10 5 g 
cm -3 ) the LPM effect becomes noticeable at k = 4000 and 
reduces cross-sections by more than an order of magnitude 
at k = 40000. However, its effect on bremsstrahlung is less 
simple. The re-calculated cross-sections in Table 4 show that 
there is an initial increase at low values of p. The explanation 
for this is as follows. For fixed values of the kinematic vari- 
ables and quantum numbers, there is a partial cancellation 
of the two terms in equation (11) which arises because the 
energy denominators are almost exactly equal in magnitude 
(« ±q) but have opposite signs. The presence of A -1 in the 
modified Green function creates a further amplitude term 
in which this cancellation is not present, so producing the 
cross-section increase. But with increasing p, values of q de- 
cline and eventually the A -1 component of the denominator 
becomes dominant. Its effect becomes marked for p > 20000 
mc. 

The presence of the LPM effect removes the simple Z 2 
scaling of cross-sections, but equations (14) and (15) allow 
us to see that A oc Z p , with -1.77 < /3 < -1.5. Thus for 



a fixed magnetic field, cross-sections become asymptotically 
dependent on Z 2+2p in the high-energy limit. 

But the possibility that zero-pressure neutron-star mat- 
ter has some degree of one-dimensional order must also be 
considered. Formation of one-dimensional chains of atoms of 
homogeneous Z, oriented parallel with the magnetic flux, is 
energetically favoured at B ~ B c (see Medin & Lai 2006) 
and such structure may be present at the surface to the 
extent that it survives the formation of the approximately 
Poisson distribution of nuclear charge produced by photodis- 
integration reactions in the reverse-electron electromagnetic 
showers. Its effect can be studied simply by reference to the 
structure factor given by equation (21) neglecting, for the 
moment, any inhomogeneity in atomic number. Consider, 
for example, the structure function for a finite linear chain 
of N a homogeneous atoms with spacing Zj — Zj-i = a s , 



An(N a qa s /2) 
sm(qa 3 /2) 



(21) 



This has maxima, S(q) = N a , when q coincides with an in- 
tegral multiple of the basic reciprocal-lattice wavenumber 
g = 2ir/a s , but has zeros at \q — g\ = g/N a . Its effect is 
primarily to significantly distort spectral shapes given by 
isolated-atom differential cross-sections. Unfortunately, it is 
not easy to be specific about changes in total cross-sections 
because they are dependent on the (unknown) structure 
function and on all other parameters, but with a change 
of variable in a particular differential cross-section to q we 
can note that the modified total cross-section is, 



(22) 



taken over the complete interval of q for the process, and 
that this differs from the S = 1 cross-section by a factor of 
order unity, not N a . 



7 CONCLUSIONS 

The cross-sections obtained here are those that, with Comp- 
ton scattering, determine electromagnetic shower develop- 
ment. For magnetic flux densities of the order of 10 12 G, 
Landau quantum numbers n ~ 10 2 are associated with 
transverse momenta p± ~ mc. Thus, on the basis of the cor- 
respondence principle, zero-field expressions are used with 
some degree of confidence. But in many neutron stars, fields 
more than two orders of magnitude larger have been in- 
ferred from the observed spin-down rate. For this reason, 
it is necessary to be aware of any field-dependence of these 
cross-sections that would qualitatively change the nature of 
shower development. In this respect, the processes that ap- 
pear not to have been considered previously are those of 
Sections 2, 4 and 5. Of these, Rutherford scattering mer- 
its little comment except to note that its cross-section de- 
creases with increasing magnetic field in the region B ~ B c . 
Pair creation is a little more complex because the ampli- 
tude (equation 8) becomes singular at the threshold values 
of k± for magnetic pair creation, beyond which it becomes 
no more than a Coulomb field- dependent correction to that 
process. Thus we are concerned with the cross-section for k± 
below the first two thresholds discussed following equation 
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(6) in Section 3. This is the explanation for the intervals of 
k± adopted for columns 3-5 of Table 3. 

Bremsstrahlung at B ~ B c is an even more complex 
process. It is of second order in the coupling of the electron 
with the electromagnetic field, as at zero field. But there 
are two significant differences. The Table 4 cross-sections 
are for an electron initially in the n = Landau ground 
state, but in the development of a shower in condensed 
matter, higher-n states are certainly populated owing to 
the very large Rutherford cross-sections listed in Table 1. 
Secondly, a pole in the amplitude (equation 11) for such 
initial states represents the possibility of real intermediate 
states of n' ^ 1. The process we refer to as bremsstrahlung 
then consists of two separate first-order processes; Ruther- 
ford scattering followed by cyclotron emission with the nat- 
ural widths given in Table 2. The initial electron momentum 
determines whether or not this is more important than the 
second-order process. An estimate of the critical value is 
easily made from the cyclotron emission rates in Table 2 
which are all of the order of 10 Mb. Although the mean 
free path Ir u for Rutherford scattering is always small com- 
pared with the radiation length ln a d which is characteristic 
of the second-order bremsstrahlung process, time dilation 
can extend the mean flight path l ce for natural decay by 
cyclotron emission to lengths greater than ln a d- Then the 
first-order processes are the dominant source of photon pro- 
duction provided the initial-state momentum is, 



V 



< 



(23) 



From the transition rates given in Table 2 and the values 
of a Rad given in Table 4, we can see that its value is not 
strongly B-dependent. Specifically, for B = B c , Z — 26, 
a Rad _ 22 bn and l Rad = 1.8 x 10~ 5 cm with inclusion 
of the LPM effect, so that for n = 4, the inequality (23) 
becomes p < 1100 mc. At lower momenta, the process of 
bremsstrahlung becomes simply a sequence of excitation 
(or de-excitation) of higher-n Landau states by Ruther- 
ford scattering and de-excitation by cyclotron emission. At 
even lower momenta, the balance between excitation and 
de-excitation changes and the typical values of the Landau 
quantum number decrease until for Ir u « l ce , that is, for 

-2- « — r u (24) 
mc c 

they are n = 1. With the cross-section 001 from Table 1, 
we find Ir u = 4.0 x 10~ 7 cm. The transition rate Tio = 
1.06 x 10~ 3 cjb gives p « 11 mc. 

The above summary of the properties of the 
bremsstrahlung and pair creation processes at high magnetic 
fields serves to show qualitatively how the characteristics of 
electromagnetic showers differ from the zero-field case. In 
the initial stages of an electron-initiated shower of about 
10 3 GeV, the Landau-Pomeranchuk-Migdal (LPM) effect is 
far more effective in reducing cross-sections owing to the 
very high density of the condensed matter at the neutron 
star surface. Therefore, shower development is displaced in- 
ward until typical electron or photon energies are reduced 
to the extent that the effect becomes unimportant. At this 
stage, electron energy loss becomes rapid, occurring within 
lengths much shorter than ln a d- Photon mean free paths are 
dependent on k±, but below the lowest threshold for mag- 



netic pair creation at k± < 2 they are a little larger than the 
value of given by the zero-field Bethe-Heitler cross-section. 

The principal interest in shower development here is in 
proton production by decay of the nuclear giant dipole res- 
onance (GDR). Cross-sections for its electro-production are 
small so that the important shower property is the distribu- 
tion of total photon track length in the neighbourhood of the 
resonant momentum, k ~ 41 mc. An estimate of this, made 
under the assumptions that bremsstrahlung occurs through 
the sequence of first-order processes decribed above and 
that the Bethe-Heitler pair creation cross-section is valid, 
was given by Jones (2010). We have found that shower de- 
velopment at B ~ B c is more complex than was assumed 
in that work, the important question being the fraction of 
photon production in the GDR region that has transverse 
momentum below the magnetic pair creation thresholds at 
k± = 2 for Landau quantum numbers n+ = n- = and 
k± = 1 + VI + 2~B for n± = with n T = 1. 

It is significant that cyclotron decay, which is dominant 
in the later stages of shower development, always leads to a 
decrease in Landau quantum number, so bringing about a 
convergence of the shower toward small k±. The only quali- 
fication here is the tendency, noted by Harding & Lai (2006) 
and confirmed in Table 2, for partial transition rates direct 
to the ground state n = to be large at B > B c . Complete 
and quantitative calculations of shower development would 
require an expansion of Tables 1, 2 and 4 to much larger 
initial-electron values of the Landau quantum number, and 
a Monte Carlo calculation using them would be a substan- 
tial undertaking which has not been attempted here. Qual- 
itatively, we can conclude that the photon track length dis- 
tribution assumed in previous work is valid at B « B c . For 
very high fields at B 3> B c , the pair creation cross-sections 
become small and the photon mean free path at transverse 
momenta below magnetic pair creation thresholds are then 
determined by Compton scattering. Therefore, we expect 
that the effect of the large cyclotron partial transition rates 
to n = will begin to reduce total photon track length in 
the GDR region of momentum at B 3> B c . 
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APPENDIX A: 
EQUATION 



SOLUTIONS OF THE DIRAC 



The functions are then, 



The solutions used here are the spinors obtained by Johnson 
& Lippmann (1949) for a free electron in a uniform magnetic 
field with the gauge A = (B x r)/2. As these authors noted, 
they can be easily expressed in terms of any complete set of 
eigenfunctions of the non-relativistic Schrodinger Hamilto- 
nian, whether in cartesian or polar coordinates. The Dirac 
equation is, 



(a ■ 7r + /3m) * = E^>, 



(Al) 



in which tt = p + eA with e > and c = 1. The solutions 
for E > and for Landau quantum number n ^ are, 



l,n 





(m + E)if!„ 

P±1pn-1 
-pipn 



*l,n- 



(m + E)ip„-i 


ptpn-1 



each multiplied by the common normalization factor G = 
2E(E + rn). The longitudinal momentum component 
is p. These solutions have eigenvalues given by equation (1) 
and are two-fold degenerate for n > 0, but they are not 
identical with the true physical states owing to radiative cor- 
rections that have been neglected. Natural line-widths, and 
level shifts arising from the electron anomalous magnetic 
moment and from radiative corrections, are all of similar 
magnitude and give physical states that are quite different 
from the above solutions (see Herold, Ruder & Wunner 1982; 
also Baring, Gonthier & Harding 2005). The distinction 
would be important if spin-polarized electron or positron 
states were being considered. However, for the cross-sections 
obtained here, which are for relativistic energies, lepton spin- 
polarization is not of interest and the fine-structure in energy 
levels can be safely neglected. In this paper, we always aver- 
age or sum over lepton spin states so that we require only the 
completeness of the Johnson-Lippmann states. These, with 
the negative energy solutions given by the transformation 
E — > —E, have been used explicitly in Sections 2-5. The 
corresponding representation of the Dirac matrices is easily 
found from equations (47)-(54) of Johnson & Lippmann and 



a- 
a 



I 

-/ 



The subscript s = ±1 labelling the spinors indicates that 
for the case n = 0, or in the limit B — > 0, they are also 
eigenstates of o z with eigenvalues as indicated. The func- 
tions i/>n are a complete orthonormal set of eigenfunctions 
of 7r 2 labelled by the Landau quantum number n. Their rel- 
ative phases are not arbitrary because they are required to 
satisfy the conditions, 



(tT x + iTTy)tp n -l = p±(ri)lp n , 

(n x - i-K y )ip„ = p±(n)ip n -i, 



(A2) 



where p± (n) = \/2nB and B is the magnetic flux density in 
units of B c . The atomic Coulomb potential has axial sym- 
metry, reflecting the shape of the electron orbitals in the 
magnetic field. (We neglect any small component with a lat- 
tice symmetry.) Thus cylindrical polar coordinates r±, <j>, z, 
with B parallel with the z ~ axis appear to be the optimum. 



ipnl 



t (i-n)/2 -«/2 



-i(l — n)4> 



2ir 



(A3) 



conveniently expressed in terms of the variable £ = r^/2r%, 
where rg = mc 2 \/B/e is the cyclotron radius. The Landau 
quantum number can be any positive integer n = 0, 1,2, .. 
and all except the n = state are two-fold degenerate. In 
addition, the quantum number I represents the infinite spa- 
tial degeneracy of all Landau states whose guiding centres 
can be at radii \J21 + lrg from the z-axis. The polynomials 
fn are, 

Jo = 1 

h = i-t 

h = /(/-l)-2i£ + £ a 
h = Z(Z-l)(Z-2)-3/(/-l)£ + 3Z£ 2 -£ 3 
U = Z(Z-l)(Z-2)(Z-3)-4Z(Z-l)(/-2)£ 
+ 6Z(Z-l)£ 2 -4^ 3 + £ 4 



(A4) 



in obvious sequence and they satisfy, 

/n(I,0 = (-O""7j(».O. 



(A5) 



which is worth noting in relation to states with I < n. 

We adopt gaussian cgs units so that the fine struc- 
ture constant is e 2 /he. Perturbations of the correct sign 
are introduced into equation (Al) by modifying the energy- 
momentum four- vector, A — > A + A and E — > E — Ze 2 V, to 
include the radiation field A and the atomic Coulomb po- 
tential for nuclear charge Z. Thus the absorption of a photon 
with wave-vector k and angular frequency u is represented 
by the matrix element of 

ea ■ A = {^f-^j ' ea • ee lk r , (A6) 

in which e is the polarization vector either parallel with, or 
perpendicular to, the vector k x B. 

The creation of a pair in states specified by 
\p±,n±,l±, s±), as in lines three and six of equation (8), 
is given by the matrix element 



/ 



r±dr±d<fidz li &l_ !n _ (p_, E- 

*-s+,n + (-p+, 



, l-)eoL ■ A 
-E+, /+), 



(A7) 



in which the positron energy is E+ > 0. Evaluation of this 
proceeds by matrix multiplication and immediately yields, 
for each set of values of s± and of photon polarization, a 
linear combination of terms each of which, apart from kine- 
matic factors, contains an integral of the form, 



/ 



I(k±, n-,l-,n+, 1+) = 
r ± dr ± d^e ik ^ cos *rp n _ :l _ (r x )^„ + , i+ (r ± ). 



These expressions, and the integrals contained in them, have 
all been evaluated numerically. Therefore, we are unable 
to give useful analytical expressions for any of the cross- 
sections calculated in this work. Evaluation is the more te- 
dious because the transverse dipole approximation is not ad- 
equate so that there is no radiative transition selection rule 
that is effective for physically significant values of k±. But 
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matrix elements of the Coulomb potential satisfy 8{n—l) — 
strictly for isolated atoms; also in condensed matter to the 
extent that components with the symmetry of the crystalline 
field are small. 

This paper has been typeset from a TjrjX/ P^rgX file prepared 
by the author. 



